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Abstract 

We show, for previously uncalculated examples containing a uniform mass in the loop, that 
it is possible to obtain complete massive one-loop gauge theory amplitudes solely from unitarity 
and known ultraviolet or infrared mass singularities. In particular, we calculate four-gluon 
scattering via massive quark loops in QCD. The contribution of a heavy quark to five-gluon 
scattering with identical helicities is also presented. 
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1 Introduction 



The computation of loop amplitudes by traditional methods is notoriously difficult involving labo- 
rious calculation. At one loop, traditional Feynman diagram techniques, augmented by computers, 
are in wide-spread use for up to four external particles. Beyond this, computations become signif- 
icantly more complicated and traditional approaches are problematic. Non-traditional approaches 
have recently been used to compute all one-loop massless five-parton amplitudes [2|, |§ and a 
variety of infinite sequences of one-loop amplitudes ||, |5|, ||, [jj . The technical developments which 
allowed for this progress include methods involving spinor helicity ||, color decompositions ||, su- 
persymmetry [ Hi, O, 1^] , string-theory (l^, [l4j , recursion relations g|, i> factorization |, [l7| 
and unitarity | lq , y, [7]. 

These developments have been mostly limited to massless amplitudes. However, in many phys- 
ical processes, such as near particle thresholds or in the production of heavy states, one cannot 
neglect masses. There is, for example, some concern that the underlying cause of the apparent 
discrepancy between the low and high energy measurements of the strong coupling constant is due 
to ignoring masses [fL9| . 

In this paper we extend unitarity techniques to obtain particular examples of massive one loop 
amplitudes directly from tree amplitudes. Unitarity has been a useful tool in quantum field theory 
since its inception. The Cutkosky rules [jl8| , |20| l allow one to obtain the imaginary (absorptive) 
parts of one-loop amplitudes directly from tree amplitudes. Traditionally, dispersion relations are 
then used to reconstruct real (dispersive) parts, up to ambiguities in rational functions. Although 
computationally more powerful than Feynman rules, the rational function ambiguity has hampered 
the use of unitarity as a method for obtaining complete amplitudes. 

Recently, it was shown that for certain 'cut-constructible' massless amplitudes - such as super- 
symmetric amplitudes - one may obtain the full amplitude solely from the cuts || J7|. This type of 
procedure has the inherent advantage of building loop amplitudes from tree amplitudes which are 
generally compact, rather than from more complicated Feynman diagrams. With this technique, 
massless one- loop supersymmetric amplitudes with an arbitrary number of external legs, but fixed 
helicity structure, have been obtained. The main new ingredient used to remove ambiguities from 



cut constructions is the complete knowledge of all integral functions |[2l] , 22 , 23 1 entering into mass- 
less one-loop amplitudes. This severely limits the analytic form of amplitudes. The cut construction 
method may be extended to full massless QCD by computing to higher order in the dimensional 



regularization parameter 24, 25]. In this way, all terms in any massless one-loop amplitude become 
associated with branch cuts which may be used to reconstruct the complete amplitude. 

In general, one is also interested in massive theories such as QCD with heavy quarks. The 
amplitudes of such theories contain logarithms that depend only on masses; such functions do 
not have cuts in any kinematic variable. This might seem to imply that one cannot obtain all 
of the terms in massive amplitudes via unitarity. However, we will explicitly show for four-gluon 
amplitudes with a massive quark loop that the well known form of their ultraviolet or infrared- 
mass singularities provides sufficient information on the coefficients of all potentially ambiguous 
functions. We also present the result for the scattering, via a massive quark loop, of five-gluons with 
identical helicities. (These amplitudes have not been computed previously.) The efficiency of the 
method for these examples may be characterized by the fact that intermediate algebraic expressions 
are about the same size as final expressions for the amplitude. Due to this feature there is no need 
for computer assistance in performing the computations presented in this paper. This may be 
contrasted with Feynman diagram calculations which frequently suffer from an exponential growth 
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in the algebra of intermediate expressions. For other five- and higher-point amplitudes, computer 
assistance is useful since the number of cuts and terms in the final amplitudes proliferates. 

In section ||, we briefly review a few of the established techniques for computing amplitudes. 
Section [3| illustrates the ambiguities that are traditionally associated with unitarity based calcula- 
tions. We then remove these ambiguities in section |I[ In sections || and [6|, we use the results of 
the previous sections to calculate all four-gluon helicity amplitudes with a heavy fermion loop. We 
also present the five-gluon amplitude with all identical helicities. Our conclusions are presented in 
section [?]. We collect various useful results in the appendices. 



2 Review of previous techniques 

We now briefly review some of the techniques which we shall use in this paper. Two techniques 
that we utilize are spinor helicity and color decompositions. The reader may consult the review 
article of Mangano and Parke [^] for further details. 

In explicit calculations with external gluons, it is usually convenient to use a spinor helicity 
basis H which rewrites all polarization vectors in terms of massless Weyl spinors \k ). In the 
formulation of Xu, Zhang and Chang the plus and minus helicity polarization vectors are expressed 
as 

where q is an arbitrary null 'reference momentum' which drops out of the final gauge-invariant 
amplitudes. The cuts of amplitudes are also gauge invariant, so one may evaluate different cuts 
using independent choices of reference momenta. We use the convenient notation 

(K\kp = (ij), (kf\kj) = [ij}. (2.2) 

These spinor products are anti-symmetric and satisfy 

(ij)[ji] = 2k i -k j . (2.3) 

This provides for an extremely compact representation of gluon amplitudes. A useful identity is 

(a\^(k;q)\b) = ±j^^(a\ [I^X^I + |^}(^|] \ h) , (2 . 4 ) 

where either (o| or \b) are spinors with four-dimensional momenta. 

For loop amplitudes, to maximize the benefit of the spinor helicity method we must choose an 



appropriate regularization scheme. In conventional dimensional regularization [27], the polarization 
vectors are taken to be (4 — 2e)-dimensional; this is incompatible with spinor helicity which assumes 
that the polarization vectors are four dimensional. An alternative is to modify the regularization 
scheme to define the polarization vectors to be four dimensional. One also defines all internal and 
external states to be four-dimensional and only continues loop momentum and phase-space integrals 



to (4 — 2e) dimensions. This is the four-dimensional- helicity (FDH) scheme [13|, which has been 
shown to be equivalent (T^j to an appropriate helicity formulation of Siegel's dimensional-reduction 
scheme |28| at one- loop. Dimensional regularization rules, which may be used in cut calculations 
of fermion loops, are reviewed in appendix |C]. The conversion between schemes has been given in 
ref. 1 12], so one may choose whichever scheme is simplest to perform calculations in. In this paper 
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we work exclusively in the FDH scheme but for fermion or scalar loops this gives the same results 
as in the 't Hooft-Veltman |2^, 27 1 scheme. 

SU(N C ) gauge theory amplitudes can be written in terms of independent color-ordered partial 
amplitudes multiplied by an associated color trace [2(| . One of the key features of the partial 
amplitudes is that the external legs have a fixed ordering. For the one-loop four-gluon amplitude 
with a fundamental representation loop the decomposition is 

Al~ loop = ffVl X>(TMi)rM2)TM3)TM4>) A 4 (a(l), a(2), a(3), a(4)) , (2 . 5 ) 

cr 

where the sum over a includes all non-cyclic permutations of the indices <r(n) and the T a are 
fundamental representation color matrices (normalized so that ti{T a T b ) = 5 ab ). We have explicitly 
extracted the coupling and a factor of fij^ from A n , where /ir is the renormalization scale. (We 
have abbreviated the dependence of A4 on the outgoing momenta kj and helicities Xj by writing the 
label j alone.) For adjoint representation loops there is an analogous decomposition with up to two 
color traces in each term. A similar decomposition in terms of 'primitive amplitudes' can also be 
performed for the case where some of the external particles are in the fundamental representation 
Q. (In refs. [|, for example, it was slightly more convenient to extract an overall factor of -/V c 
from the color ordered amplitudes, A n , but since we only deal with fundamental representation 
loops in this paper there is no need for this.) 

An essential feature of this decomposition is that the partial amplitudes, A n , are described by 
diagrams with a fixed cyclic ordering of external legs. In this paper we deal exclusively with such 
amplitudes. 



3 Naive application of Cutkosky rules 

In this section we demonstrate how a naive application of the Cutkosky rules can miss contributions 
to a complete one-loop amplitude. Traditionally, one first computes the imaginary^ (absorptive) 
part of an amplitude by evaluating a four-dimensional phase-space integral. (The phase-space inte- 
grals are not ultraviolet divergent and therefore converge in four-dimensions.) The real (dispersive) 
part is then reconstructed via a dispersion relation. We perform the cut calculation a bit differ- 
ently; since we want the complete amplitude, it is convenient to replace the phase-space integral 
with an unrestricted loop momentum integral which has the correct branch cuts (f| [7| . This allows 
us to simultaneously construct the real and imaginary parts. The ambiguities that we encounter 
are equivalent to those of the more conventional approach. 

We will compute cuts of massive amplitudes in all possible channels in an attempt to reconstruct 
their functional form. We consider the amplitude, not in a physical kinematic configuration, but in 
a region where exactly one of the momentum invariants is taken to be above threshold, and the rest 
are negative (space- like). In this way we isolate cuts in a single momentum channel. Furthermore, 
following refs. || f?j, we apply the Cutkosky rules at the amplitude level, rather than at the diagram 
level. This is advantageous because amplitudes are generally much more compact than diagrams. 

Consider the s-channel (s = (k\ + /C2) 2 ) cut of the four-gluon amplitude pictorially represented 
in fig. |l]a, and given by the phase-space integral 

3 By imaginary we mean the discontinuities across branch cuts. In Feynman diagrams these are unfortunately the 
real parts. 



3 




Figure 1: The s- and i-channel cuts of a one-loop four-gluon amplitude. The cut lines represent 
either a fermion or scalar. 



i Disc Ai(l, 2, 3,4) = / -f^-j 2tt5 (+ ^ (i\ - m 2 ) 2n6^ {l\ - m 2 

s-cut J (27r) 4 



(2tt) 4 v 1 ' w 7 (3.1) 

x AT e (~ii, 1,2, £ 2 ) 4 ree (-^2, 3, 4,4) , 

where £i = p and £2 = P — k\ — k 2 , m is the mass of the particle in the loop, <5( + ) is the positive 
energy branch of the delta-function and 'Disc' means the discontinuity across the branch cut. Note 
that, in performing the sewing of the tree amplitudes, we have maintained the clockwise ordering 
of the legs, as required by the color decomposition. 

Since we wish to reconstruct the full amplitude and not just the imaginary (or absorptive) part, 
it is convenient to replace the phase-space integral with the cut of an unrestricted loop momentum 
integral ||. This is accomplished by replacing the <5-functions, which impose on-shellness, with 



propagators: we replace eq. (3.1) with 



^4(1,2,3,4) = f 7^4 ^ir(-4,l,2,l 2 )^4 M (-43 1 V 1 ) (3.2) 
s-cut J (27r) 4 i\ — m z 12 — m z s— cut 

As indicated this equation is valid only on the s-cut. A similar equation holds on the i-cut (where 
t = (k 2 + ^3) 2 )) depicted in fig. [Tjh>. Up to potential ambiguities in rational functions, one obtains 
the full amplitude by combining the two cuts into a single function. 

For example, consider the amplitude A^(l + , 2 + , 3 + , 4 + ) with a massive scalar^] in the loop. This 
amplitude is useful for illustrative purposes because it is particularly simple. The tree amplitudes 

4 In this paper use the the conventional normalization that the number of physical states in a fundamental repre- 
sentation complex scalar is two. In refs.ji|, [?], ^[ [| a supersymmetry preserving normalization with four states in a 
complex scalar was used. 
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entering the two sides are easily computed from Feynman diagrams and are given by 

[12] 



AT e (-£ u l + ,2 + ,e 2 ) = im 



4 ree (-^2,3 + ,4+,£ 1 ) = im 



(12) {(h - h) 2 - m 2 ) ' 
2 [34] 
(34)((^ 2 -A:3) 2 -m 2 )' 



(3.3) 



where legs t\ and £ 2 represent the cut scalar lines. From eq. ( |3.2|) the scalar loop contribution to 
the four-gluon amplitude is 



,4 4 (l + ,2 + ,3 + ^ 



m 



-cut 

[12] [34] f d 4 p 



1 



(1 2) (34) J (2vr) 4 (p 2 - m 2 ){{p - k{) 2 - m 2 ){{p -k x - k 2 ) 2 - m 2 )((p + k 4 ) 2 - m 2 ) 
, m 4 [12] [34] D=4 
'16^ 2 (12) (3 4) 4 .-cut' 



s— cut 



(3-4) 

where I®~ 4 is the scalar integral box function givenF] in eq. ( A. 6 ) of appendix [A]. The function 
J 4 D=4 contains branch cuts since it contains logarithms and dilogarithms. At this stage the equality 
of both sides of the equation is only for the cut in the s-channel. 

The t-channel cut is similar and may be obtained by relabeling the legs by 1 «-> 3. Using the 
identity 

[2 3] [41] [12] [3 4] 



(2 3) (41) (12) (3 4) 
the full amplitude is found to have the form 

im 4 [12] [34] tD=a 



A 4 (l + ,2 + ,3 + ,4+) 



16vr 2 (12) (34) 



I 4 + rational , 



(3.5) 



(3.6) 



since it has the correct branch cuts. The rational terms have no cuts. Indeed, if we compare this 
to the correct result for the amplitude 



.4 4 (l+,2+,3+,4+) = -L- MlL^ =4 _ l\ 

12) 34) V 4 6/ 



(3.7) 



16vr 2 (12) (34) 

we explicitly see that the last term was missed by our naive application of the Cutkosky rules. We 
note that setting m = in eq. (|3.7|) gives the correct massless scalar loop result [13j. 



4 Fixing the ambiguity 

We now demonstrate that the last term in eq. ( |3.7[ ), which is cut-free, may be obtained from the 
first one with little additional work but with a more careful application of unitarity. 

4.1 Cuts with dimensional regularization 

The Cutkosky rules are valid in arbitrary dimensions — that is they may be used with dimensional 
regularization. By working to all orders in the dimensional regularization parameter, e = (4— D)/2, 

5 The integral given in eq. (|A.6|) is actually dimensionally regularized, but there is no difference through O(e ). 
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the number of potential ambiguities is greatly reduced. In a massless theory, it is clear that there 
are no polynomial ambiguities in dimensionally regularized expressions. All terms necessarily must 
contain cuts since the dimensions are shifted by — 2e and must therefore be proportional to factors 
of the form (— K 2 )~ e , where K 2 is some kinematic variable. For example, in the massless case, by 
working to higher order in e, the 1/6 in eq. ( [3.7| ) appears as 

which contains cuts at 0(e) linked to the 1/6. Indeed, this observation may be used to compute 



the rational part of massless amplitudes directly from the cuts [24, 25 1. 

For massive theories this simple argument requires modification, because of the presence of addi- 
tional dimensionful parameters (the masses). A factor of m~ 2€ can in principle soak up dimensions 
to generate terms without cuts. We will return to this issue at the end of this section. 

As a first step in removing the rational function ambiguities, we modify the calculation of 
the previous section to be valid to all orders in the dimensional regularization parameter, e; as 
we shall see this automatically removes most ambiguities. We consider again the cut integral in 
eq. ( |3.4| ) and fig. |l|a, except this time we replace the four-dimensional loop momenta, l\ and £2, 
with (4— 2e)-dimensional ones. A simple way to implement this is to modify the on-shell condition 
on the cuts to l\ — li 2 = m 2 and £\ — li 2 = m 2 , where t\ and I2 are left in four-dimensions, but Li a 



is a vector representing the (— 2e)-dimensional part of the loop momentum. (See Appendices A. 2 
and ^|.) With this notation, it is easy to see that the dimensionally regularized forms of the tree 
amplitudes are obtained from eq. (3A) by shifting m 2 — ► m 2 + li 2 . 



Thus, the dimensionally regularized expression for the s-channel cut replacing eq. ( |3.4| ) is 

[12] [3 4] f d 4 p d~ 2e Li (m 2 + Li 2 ) 2 



a scalar /i + q+ q+ a+\ 

[12] [3 4] f 

Ai (1 ' 2 ' 3 ' 4 K -cut" (12) (34)7 



-cut (12) (3 4) J (2vr) 4 (2?r)- 2e V 4 



(4.2) 



-cut 



where we have explicitly separated out the integration over the (— 2e)-dimensions. The box denom- 
inator is 

£> 4 = (p 2 ~ ii 2 ~ m 2 ){(p - ki) 2 -ix 2 - m 2 )((p -k x - k 2 ) 2 - li 2 - m 2 ){(p + k 4 ) 2 - li 2 - m 2 ) . (4.3) 

In obtaining this form we used the standard prescription that a four-dimensional vector is orthog- 
onal to a (— 2e)-dimensional one. This representation of dimensional regularization follows the one 



used in its original derivation [29] and has recently been used by Mahlon in his recursive approach 
151. 



As discussed in appendix [A.2| , integrals with powers of li in the numerator may be expressed in 
terms of higher dimension loop integrals, 

2V1 = „'/ 1 ( a \2—€ f di P d if 12 )" _ _ Afo _ s . . . / _ -, _ n, T D=2r+A-2e 



(4.4) 

where I®~ 2r+i 2e is the scalar loop n-point integral function evaluated in (2r + 4 — 2e)-dimensions 
and V n is the appropriate product of denominators. Thus, the correct s-channel cut to all orders 
in e is 

A scalar /-i + o+ 9+ A + \ * [1^] [34] 

A * (1 ' 2 ' 3 A ) s-cut = (4^-Ml2)(34) g4 ,-cut' ^ 
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where (see appendix A,4| ) 

K 4 = / 4 [( m 2 + M 2 )2] = m 4jD=4-2e _ ^ jD=6-2e _ ^ _ fi) ^=8-2* 

Once again the t-channel cut may be obtained by relabeling k\ <-»• &3, so that the function 

2+3+ 4+) = (4 ; )2 _ £ g| || ^4 + f(m 2 ) , (4.7) 

has the correct cuts in both the s and t channels. The function f(m 2 ) has no branch cuts in any 
kinematic channel; it has a form ~m~ 2e , but may contain rational functions of kinematics variables. 

4.2 Fixing the remaining ambiguity 

In appendix [B] we show that, for the four-gluon amplitudes considered in this paper, the following 
procedure uniquely fixes all cut-free functions. 

1. Add in the tadpole function I\ with a coefficient adjusted so that the quadratic ultraviolet 
divergence vanishes. (The quadratic divergences cancel in gluon amplitudes as discussed, for 
example, in ref. |20||.) This may be implemented by dropping all explicit 1\ functions, then 



substituting J2 — ► J2 + I\ and J 2 6 2e — > 1% 6 2e — h ■ These integral functions are defined 



in eqs. (Q, ( gig ) and flCT|) 



2. Add in the bubble function ^(0) with a coefficient adjusted so that the ultraviolet divergence 
of the amplitude agrees with the well known one. The ultraviolet divergences are 

= ^ 2 r(1 + 6) (n-2)(H^£_!V_^) Ar (48) 

n UV singular (4tt) 2 ~" £ V 'V 6 3 12/ " V 7 

for a theory with A^ c colors, fermion flavors and n s complex scalars. The first contribution 
in the parenthesis is from a gluon loop and the remaining two from fermion and scalar loops. 

We will use these rules for fixing all cut-free functions in the amplitudes presented in this paper. 

As an alternative to the second step, we may use the known infrared behavior as the masses 
vanish. For massless QCD, the dimensionally regularized infrared divergences have been categorized 



in refs. [3C]. This may be used to fix the coefficients of all ln(mj) terms as mi — > 0. In this limit 



the matching condition for infrared singularities is 

ln(mf ) ->- + ... (4.9) 

In general, the infrared properties are more powerful since we can use them to fix coefficients when 
there are multiple masses present in an amplitude. 

For the amplitudes considered here, we find that the second step is equivalent to adding in the 
cut-free self-energy diagrams on external legs. For each external gluon leg the bubble contribution 
is A tTee F/2, where 
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and 

= - ^y-* n \e e) m ~ 2e = -I Wr* l2(0) • 

corresponding to scalar and fermion loop contributions. 

For the A| calar (l+, 2+, 3+, 4+), the two steps are trivial because ^4 tree vanishes and there are no 
ultraviolet divergences. This means that f(m 2 ) = m eq. ( gjD so that 

4 scalar /•■!+ o+ q+ 4+^ _ j t 1 2 ] [ 3 4 J y , (a -iq\ 

is the complete amplitude valid to all orders in e. 

The last two terms in defined in eq. Q4.G| ), are proportional to e and will contribute through 
O(e ) only if the integrals contain an e _1 . The integral I^~ G ~ 2e is finite, but I®~ 8 ~ 2e has an 
ultraviolet divergence. Extracting this divergence either by direct loop integration, or by use of the 
integral recursion relation ( |A.12| ), we have 

-6(1 -6)/f= 8 - 2£ = -^ + 0(e), (4.13) 

which gives the missing rational term in eq. (U). Hence, eq. flOg) agrees with the correct result 
( |3.7| ) and contains all rational functions through O(e ). Since the ultraviolet divergences of the 
integral functions are independent of masses, this rational contribution is identical to the contri- 
bution in the massless case, through 0(e). This feature holds for all the amplitudes presented in 
this paper. 



5 Fermion loop amplitudes 

In this section we consider fermion loop amplitudes. Fermions are slightly more complicated than 
scalars as one must deal with Dirac algebra, a summary of which is presented in appendix |C^. 
Apart from this, the techniques used in the previous section for the scalar loop carry over to 
fermion calculations. 

5.1 The A 4 (l + ,2+,3 + ,4+) amplitude 

As a first example, we reconsider the all plus helicity amplitude but with a massive fermion loop 
instead of a massive scalar loop. In this case the tree amplitudes appearing on both sides of the 
s-channel cut (see fig. ||a) are 

Ar%-Lll + ,2+,Lf)=iW (-L^ + m)\-L 2 ) 



(12) (L.-k^-m 2 ' 

treef / + + f) [34] (-LakfQt + rrQI-L!) 
A 4 (-L 2 ,3 ,4 ,^i)-« (34) {L2 _ h)2 _ m 2 > 

where u>± = ^(1 ± 75), capitalized momenta are (4 — 2e)-dimensional and lower case momenta are 
four-dimensional. (See appendix [A|.) In these tree amplitudes we suppress the distinction between 
fermion and anti-fermion spinors: u(p) versus v(p). We find this convenient since it provides a 
uniform notation for both types of spinors. When the momentum follows the direction of the 
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fermion arrow it represents a fermion state. Conversely, when the momentum flows against the 
arrow we have an anii-fermion state. However, using this convention an overall sign error is made 
for the contribution of a vv pair; u(—p)u(—p) = —p + m = — J2 v (p)v(p)- Since this is an overall 
sign, we may easily correct for it by hand. Alternatively, one may work with both u and v to obtain 
identical results. 

The method for obtaining the cuts is identical to that of the scalar case in eq. ( |4.2D , except that 
instead of (fj? + m 2 ) 2 the numerator is initially 

-(-Li\u}+(d + m)\ -L 2 )(-L 2 \u + (d + m)\ -L x ) = 

(5.2) 

- tr[(fa + i/i - m)uj + (^ + m)(f 2 + ^ - m)uj + (^ + m)J , 
where the explicit overall minus sign compensates for our suppression of the distinction between 



fermions and anti-fermions. Since ^ (like m) commutes freely (see eq. ( C.7 )) with the helicity 



projection operator, cj+, we find that the fa in the above trace cancel. For example the term, 

• • • oj + {i/l + m)fa oj + ■ ■■ = •■■ + m) u + fa u + ■ ■ ■ = . (5.3) 

The trace thus reduces to a combination of ^ and m containing terms. We readily evaluate what 
remains as, 

tr[(^ - m)uj + {i/L + m){i/i- m)uj + (ft + m)] = tr[w + (^ 2 - m 2 ) 2 ] = 2(/i 2 + m 2 ) 2 , (5.4) 

so that 

,fermion n + 9 + o+ , + N _ 2^ [1 2] [3 4] 

A 4 (1 ,2 ,3 ,4 )- (47r) 2^ (12)(34) K 4, (5.5) 

where is defined in e q. (f4.6| ). We have fixed the cut-free function, by following the same steps as 
for the scalar loop case (|4,12j ). Thus for this helicity, the contribution for a Dirac fermion is minus 
twice that of a complex scalar given in eq. ( 4.12[ ). This is the correct result as may be verified from 



a supersymmetry identity [10|. In summary, for this helicity, we have obtained the correct result 



for a massive quark loop, with little more work than with a naive application of unitarity. 
5.2 The A 4 (l ,2 ,3+,4+) amplitude 

Having calculated the simplest of the closed fermion loop sub-amplitudes consider now one of 
the more difficult ones. The ^(l - , 2~, 3 + , 4 + ) sub-amplitude does not have an s <-> t exchange 
symmetry. Accordingly, we must explicitly calculate both cuts and combine the information to 
obtain the complete result. Also, this amplitude contains non-vanishing cut-free functions whose 
coefficients will be fixed from the ultraviolet divergences. 

The easier cut is in the s-channel, so we calculate it first. The required tree- level amplitudes 

are, 

.tree, rf n rA _ A 1 2 ) W^llM 

4 ( 1,1 ' ' L2) - Z [12] (Li — ki) 2 — m 2 ' 

t / + + f [34] {-L 2 \^ + m)\-L l ) 

A 4 (-L 2 ,3 ,4 ,^iJ-* (34) {L2 _ h) 2_ m 2 ■ 
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As before, these combine to give a contribution to the s-cut, 



-cut 



s 



drp d 2e p 
(2tt) 4 (2vr)- 2e 

{-Lx\u-(iA + m)\ -L 2 ){-L 2 \u + (d + m)\ -L x ) 



where t\ = p, t 2 = p — k\ — k 2 and we use (see ref . |26l] ) 



V, 



( 12 ) [34] _ J (IT = ■ t t + + 

[12] (34) s (12) (23) (34) (41) s 4 1 ' ' ' J ' 

Concentrating on the numerator Dirac algebra under the integral, we create a trace, 

(-Li|w_(^ + m)| -L 2 )(-L 2 |w+(^ + m)\ -Lx) = 

tr[(/x m)uj-{ifL + m)(i 2 + ft - tn)uj + (ft + m)) . 



s— cut 

(5.7) 



(5.8) 



(5.9) 



Unlike the all plus case, the helicity projection operators do not cancel the four-dimensional mo- 
menta. Instead they cancel the (d — m) contributions from the spinors on the cut. The above trace 
reduces to, 

tv[ujx(l^ + ™%{ft + m)) = {p 2 + m 2 )ti[LuJJ 2 } = (p 2 + m 2 ) 2£i • £ 2 

0/ 2 1 2\2 f 2 1 2\ (O.lUj 



where we have used the fact that on the cut t\ 



t 2 = p 2 + m 2 . This integrand exhibits a 



supersymmetry decomposition as the first term is identical (up to an overall factor of —2) to the 
contribution of complex scalar loop. 

Performing the integration, the s-cut of the amplitude becomes 



^fermion 



(1-2-3+4+) 



-cut 



(4vr 



1 A ^( 2 -K 4 + tJ 4 

V s 



s— cut 



(5.11) 



where J4 and K4 are defined in eqs. (A. 24) and (A. 25). 

The i-cut, depicted in fig. |]b, is more complicated. It may be constructed from the following 
two tree amplitudes, 



A t j ee (-L{,2-,3+,L f 2 )=ie 2 --£ 1 
A i i ee (-L(,4+,l',L f 1 ) = ie+-e 2 



(-Lx\tl 



?l — k^ 2 — p 2 — m 2 
(-L 2 \fc\-Lx) 



(5.12) 



(h - h) 2 - p 2 



where the reference momenta qi for the polarization vectors £j are are q 2 = k$, = k 2 , q<± = k x , 
and qi = k^. This is a particularly good choice of reference momenta as each tree amplitude is 
described by a single Feynman diagram. As depicted in fig. |l|b, for the i-channel the cut momenta 
are 

ix=p-k x , £ 2 =p + k 4 . (5.13) 

This representation of the tree amplitudes may be directly obtained from Feynman diagrams after 
using the spinor helicity form of the polarization vectors. (As discussed in appendix |C], identity 
Q2.4|) may be used on the j, only when it is adjacent to a four-dimensional quantity.) 
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Using eq. Q3.1| ) and combining the two tree amplitudes into a trace, the i-cut becomes, 

(2tt) 4 (27r)" 2e 

e 2 • £i£4 • ^2tr[(-/x - i/l +m)^{-f. 2 -^+m)^J 



4 ermion (l-,2-,3 + ,4+)| i t = -/ 



Pi 



t— cut 

(5.14) 



(5.15) 



The numerator may be rewritten as 

N = -e 2 ■ £184 ■ ^ 2 tr[(/i +//L -m)^ 3 (/ 2 +^ -m)^] 
= -2A/i - M 2 , 
where 

A/i = 4e 2 -£i£4 -^£3 • 4ei ■ 4, Af2 = 2te2-£ie A -e 2 ei ■ e 3 . (5.16) 

We have dropped terms with an odd power of /i a (see appendix [A] for the relevant discussion). J\f\ 
may be recognized as the numerator of a cut complex scalar loop, while the second is the additional 
piece for a fermion loop. Thus, this integrand again obeys a supersymmetry decomposition. 



To evaluate the amplitudes, we use spinor helicity. From eq. (2.1), and the choice of reference 
momenta given below eq. ( |5 . 1 2| ) , the polarization vectors satisfy 

, (12) [34] (4 + |A|l + ) , (3+|/i|2+) 

£ 1 ' e 3 — o\ r-. .1 > e l - *2 



(23) [14]' * ^2[14] ' ' v^[23] ' 

63 ^2(23) ' " 4 ' £2_ ^2(14) • 

It is convenient to evaluate separately the two numerator terms of eq. ( |5.15 ). The first corresponds, 
up to a factor of —2, to the contribution of a complex scalar loop, whose numerator is 

= (4+|A|l+) (3+1^12+) (2-|/ 1 |3-)(l'|/ 1 |4-) 

[14] [2 3] (2 3) (14) (518) 

= i(2-|/ 1 |3-) 2 (4+|/ 1 |l+) 2 . 

By multiplying and dividing by (34) 2 [1 2] 2 , J\f\ can be rewritten in terms of a trace as 

M = p<2~|/i |3-) 2 ||J<4 + |/ a |1 + ) 2 |^ = ^< ee (l-,2-,3 + ,4+)tr + [^34a2£34a2] , 

(5.19) 

where we have defined l\=£, and for brevity represent the external momenta by their indices and 
neglect slashing the momenta inside the trace. (We use the notation, tr±(. . .) = tr(cj± ■••)•) Since 
tensor box integrals are generally difficult to evaluate (without the aid of a computer), we will 
avoid such terms. To simplify the trace, we commute the is within the trace towards each other. 
This generates more simple integrals, with inverse propagators or fewer powers of loop momentum 
in the numerator: 

tr+[e34£12£34£12] = tr+ [£34£12£342] (11) - tr+ [£34£12^341] (£2) + tr+ [234^12^341] [pi 2 + m 2 ) 

= tr+[£34£l2&42] (£ 1) - tr + [£34£12^341] (£2) + tr+ [2342^341] (fi 2 + m 2 ) (£1) 
- tr + [2341£341] {^ 2 + m 2 ){£2) + tr+ [23412341] (/? + m 2 ) 2 . 

(5.20) 
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Using 



(£l) = 2£-k 1 = p 2 -^ 2 -m 2 , (£2) = 2e-k 2 = -(p-k 1 -k 2 ) 2 + lJ- 2 + m 2 , (5.21) 

this becomes 

tv + [£3All2£342](p 2 - /J 2 - m 2 ) + tr+[^34^12£341]((p - h - k 2 ) 2 - fi 2 - m 2 ) 
+ tr+ [2342^341] (/i 2 + m 2 )(p 2 - fj, 2 - m 2 ) + tr+ [234 W341] (/i 2 + m 2 )((p -k x - k 2 ) 2 - /i 2 - m 2 ) 
+ tr+ [23412341] (/i 2 + m 2 ) 2 , 

(5.22) 

which exhibits inverse propagators in the numerator. After cancelling these inverse propagators, 
against those present in eq. ( 5.14 ), we obtain tensor triangle and scalar box integrals. 

One way to evaluate the tensor integrals is via Feynman parameterization. A convenient way 
to do this is to have a uniform labeling of Feynman parameters for all diagrams: for the triangles, 
the Feynman parameter corresponding to the missing internal propagator is set to zero. With this 
convention, the Feynman parameter shift for all diagrams is, 

£ = p — k\ = q— k\a\ + /c 2 (a 3 + a 4 ) + fc 3 a 4 , (5.23) 

where q is the shifted loop momentum. For example, the first term in eq. (5.22) corresponds to the 
triangle integral 

J_ .tree f d A P d~ 2e \i tr+ [1341121342] 

a 5 * 4 J (2vr) 4 (2vr)- 2e ((p - h) 2 - /i 2 - m 2 )((p - k x - k 2 ) 2 - /i 2 - m 2 )((p + k 4 ) 2 - /j 2 - m 2 ) ' 

(5.24) 

Performing the Feynman parameter shift in eq. ( PT23D , but with ai = since the corresponding 
propagator cancels, and integrating over q and \x we obtain (see appendix |~~ 

2i ,t,«, f l * „, , t dt 1 d ~ 2 'f> tr + ( g 34(<, + /i 2 03)125342] 



(2vr) 4 (2vr)- 2e (q 2 - a 2 a 4 s 23 - m 2 - /j 2 ) 3 

-(^< ee ^tr + [ 7 ,342 7 ^342] S 4 1 )' D = 6 - 2e [a 3 ] ^ 



= 0, 

where we use = 0, momentum conservation and 

Iv^m^nl = 4Ato " k n i ^/u^m^n^ll = ~'^F'lF'n£m • (5.26) 

(Note that the Lorentz indices of the 7-matrices are four-dimensional. This means that the required 
7- matrix identities are the standard four-dimensional ones.) As defined in appendix A3, the 
integral function I n [ a i] is a Feynman parameterized form of the loop integral with a factor of aj in 
the numerator. The superscript on I'tp indicates the triangle obtained from the box by removing 
the propagator between legs 4 and 1, as discussed in appendix |A]. 

Similarly, evaluating all the terms in eq. ( |5.22| ) we obtain 



^r ar (i-,2-,3 + ,4 + ) = — ^4 ree 4(-i tr +[ 234 > 12 ^ 341 ]4 3),D " 6 " 2e [«4] 

+ tr+ [23412341] } [a 4 ] - tr+ [23412341] £T 4 ) , 
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where jf } [a 4 ] = m 2 ^ [a 4 ] - eif } ' D - 6 " 2e [a 4 ]. (cf. eq. 

Using momentum conservation to evaluate the traces we have 

Af* = ^4 ree (4 3),D=6 - 2£ N + ^ 3 (3) K] - ^4) | t _ cut ■ (5.28) 

The 75 in tr+ vanishes for a four-point function, because there are only three independent momenta 
to contract into the Levi-Civita tensor. 

As a final step, we may convert the single parameter triangles to scalar integrals using the 
integral relation ( A.22j ), which explicitly yields 

if [a4 ] = \ ip) , if' D=6 " 2 > 4 ] = \ ja.3)^=6-2e > (5>29) 

where we have dropped those bubbles without a t-channel cut. This gives the final form of the 
t-cut as 



^scalar = _}_ ^tree / 1 J (l,3) )0 =6-2e 1^1,3) _ t N (5 _ 3Q) 



Combining the t- and s-cuts in eq. ( 5.30| ) and minus half of the first term in eq. ( |5.11 ), we find 
that the full scalar loop amplitude is of the form, 

^ calar (l-,2-,3+,4+) = j^Ar e (]l^ D=6 - 2t + - V) +f(m 2 ), (5.31) 

where f(m 2 ) is function of the masses which does not have a cut. 



We still have to evaluate the uniquely fermionic term in eq. (5.15), A/2, which combines with 
the scalar loop contribution to form the complete fermion loop. Since it contains only two powers 
of loop momentum in the numerator, instead of four, it is simpler to evaluate than the above scalar 
loop term. The result of such an evaluation is 



4 ermion (l- 2~ 3+ 4+) = -2Af alar (l~,2~,3+4 + ) _ _L_ A ^ ( t j 4 _ 1^) + g ( m 2 ) . (5.32) 

(471") ' 



Finally, we fix the form of f(m 2 ) and g(m 2 ) by correcting the ultraviolet divergences with the 
two-step procedure in section §.2 to obtain the final result 

Ar lar (l-,2-,3+,4+) = ^_^ee(l (/ D=6-2e W _ ^ + I(J 2 ( t ) _ t Ki - I/ 2 ( )) , 



(47 

1 

(4t0 



4 ermi ° n (l~>2-,3 + ,4+) = -2vl| calar (l-,2-,3 + ,4+) - j—^-A^ (t J 4 - i 2 (i) + 2i 2 (0)) , 

(5.33) 

where A^ ee is defined in eq. (5.8). These expressions for the amplitudes are valid to all orders in 
the dimensional regularization parameter, e. 

These amplitudes are bare ones before subtraction of the ultraviolet singularities. To obtain 
the MS renormalized amplitudes (in the FDH regularization scheme), we subtract the ultraviolet 
divergences 

- — — — r ^,"t e ^ ree - , for fermion loop, 

3 (4?r) 6 (5.34) 

- — ' ^ r ^2-1 \ ' fOT SCalar lo °P' 
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where n = 4 is the number of external legs. 

After appropriately redefining the coupling constant [12], the results in eq. ( 5.33 ), are identical 
in the 't Hooft-Veltman scheme (where observed particles are in four-dimensions, but unobserved 
ones are continued to D = 4 — 2e). (The gluon loop contribution is, however, shifted |p^3| .) One may 
also convert to the conventional dimensional regularization scheme, by accounting for the external 
e-helicities 



Whilst unitarity unambiguously fixes the correct sign for the cut of any loop amplitude, it is 
frequently convenient to ignore the overall sign of a cut as it may be determined in the final stages 
of a calculation. By comparing integral functions with cuts in multiple channels we can usually 
fix the relative signs of the separately calculated contributions. The overall sign of the amplitude 
can be obtained by comparing with known results in the massless limit, or in the case of five- or 
higher-point amplitudes, factorization properties are a sufficient constraint. Alternatively, these 
consistency requirements may be used as checks. 



6 Remaining amplitudes 

The remaining four-gluon helicity amplitudes are straightforward to calculate using the above 
techniques. For completeness we quote the results: 

^scalar/-!- o+ o+ 4 +\ [ 

4 { ' ' ' ; (47r) 2 -^ [12] (23) (34) [41] u 

x ( J 3 (s) + J 3 (i) - ^(J 2 (j) + h) ~ S -^(Mt) + h) + £ J 4 + K 4 ) , 

V su tu s z t z lu J 

A scalar (1 - )2+)3 - )4+) = (^2- J 3 (i) + ^^J 3 ( S ) - ^/ a (t) - £l 2 (s) 



(47r) 2 e V U 3 ' U 3 M 3 M 3 

f (I?= 6 ~ 2 *(t) - h) - ±(I?=^(s) - h) - 4(J 2 (t) + I X ) - - 2 {J 2 (s) + h) 
tu su u z u z 

u 2 3 



(6.1) 



where u = —s — t and 



The remaining fermion loop amplitudes are 

A 4 ermion (l-,2+,3 + ,4 + ) = -2^ calar (l",2 + ,3 + ,4+), 
A 4 ermion (l-,2 + ,3-,4+) = -2^ calar (l',2+,3-,4+) 

- jA^ A T e (-h{s) + S - h{t ) - ^ J 4 - *I? = 6 "^ + 2/ 2 (0)) . 
(47rj z e \u u u u > 

(6.3) 

These amplitudes explicitly exhibit the supersymmetry decomposition discussed in refs. [|], |TT]. [|, 
since the fermion loop amplitudes contain pieces which are naturally expressed in terms of the 
scalar loop amplitudes. 

We have explicitly checked that the four-point amplitudes obtained via unitarity are in full 
agreement with a more conventional diagrammatic calculation. (The forms obtained via a conven- 
tional Feynman diagram calculation are in general somewhat more complicated than those obtained 
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via unitarity; the two forms may be compared using the integral recursion formulas ( A.12p ). As 



before, to obtain the renormalized amplitudes from the above bare amplitudes we subtract from 
the amplitudes the quantities in eq. J5.34 ). 



The methods we have presented here can also be applied to computing higher point functions. 
For example we have obtained the simplest of the five-gluon amplitudes: 

4 ermion (l + ,2+,3 + ,4 + ,5 + ) = -2^ s 5 calar (l + ,2+,3 + ,4 + ,5+) 
2i 1 



(12) (2 3) (3 4) (4 5) (51) (4^- 

— - [s 23 S34tr_ [5124]^ 1} + s 34 S45tr_ [1235]#f } + s A5 s 5l tr„ [2341]k| 3) 



tr 5 [1234] 

+ S 5 iSi 2 tr_[3452]i4 4) + Sl 2 S23tr_ [4513]^ 5) - Si 2 S23S34S45S51^5 

where Sij = {ki + kj) 2 , tr$[A] = tr[75^4] and tr_[A] = ^tr[(l — 75) A], The function K± is obtained 
from eqs. ( |A.S| ), ( |A.12[ ) and ( [A.25| ). We may use eq. ( A.9| ) to express K$ in terms of K± . The (i) 



label on K A follows the same convention as that used for I±. Extending the discussion in appendix [B| 
to this case, it remains true that the only cut free functions that can appear are I\ and I 2 (0). Since 
the corresponding tree amplitude vanishes, this amplitude must be infrared and ultraviolet finite; 
the coefficients of these two cut-free functions vanish as was the case for eq. flOp . 

A few simple consistency checks may be performed on this amplitude: in the massless limit 
through O(e ) it reproduces previously calculated results Ej. Furthermore, it is not difficult to 
check that it satisfies appropriate limits as any two external momenta become collinear g|, g 0. 
(In performing the check to all orders in e, integral identities of the form discussed in section 5.2 
of ref. |l7j are required.) 

In obtaining these results, no computerization was necessary. For other helicity structures, or 
different external particles, the amplitudes become sufficiently complicated that computer assis- 
tance is useful. 



7 Summary and discussion 

In general, it is far easier to build amplitudes from their analytic properties than from diagrams. In 
order to use unitarity (l8|, ^] as a computational tool for complete amplitudes, one must maintain 
control over all potential cut-free contributions. A powerful constraint is that the functions entering 
into an amplitude are not arbitrary, but must be composed of a definite set of integral functions. 
This has recently been used to obtain certain infinite sequences of one-loop amplitudes || . With 
unitarity, the basic building blocks of one-loop amplitudes are tree- level amplitudes, which generally 
are far more compact than individual diagrams, especially in a helicity representation. Unitarity is 
also esthetically appealing in providing a way to build further amplitudes using on-shell quantities. 

In this paper, we extended the unitarity method of refs. H f?j to obtain particular examples of 
one- loop amplitudes with a uniform mass in the loop. We have computed all four-gluon helicity 
amplitudes with massive fermion or scalar loops. These amplitudes have not previously been 
computed, and are useful for investigating the heavy quark thresholds in elastic glue-glue scattering. 
If unexpected structure were to be observed in the experimental jet cross-sections at heavy quark 
thresholds, a detailed analysis of this effect would be warranted. 
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We have proven that all potential cut-free functions that may appear in these amplitudes can 
be fixed with a knowledge of the ultraviolet or infrared divergences. By performing the calculation 
to all orders in the dimensional regularization parameter, only a few cut-free functions can appear 
in the amplitudes. These cut-free functions are integrals with no dependence on external kinematic 
variables. For the examples treated in this paper with a uniform mass in the loop, computing to all 
orders in e (where D = 4 — 2e) takes little extra work as compared to working through 0{e ). The 
reason is that the e components of the loop momentum act as a uniform mass in the loop that is 
then integrated over ||] . We have also demonstrated that the coefficients of the cut-free integral 
functions may be fixed using well known ultraviolet or infrared divergences of gauge theories as 
'boundary conditions'. 

The amplitudes obtained in this paper had a number of generic features. Firstly, supersymmetry 
relations [ll]] between contributions to the amplitudes were identifiable at the level of integrands. 
This may be contrasted to conventional Feynman diagrams where no such simple identification is 
possible. This allowed us to compute the scalar loop contribution as a byproduct of the fermion 
result. A second property was that rational functions (such as in eq. (|4.13 )) appearing in the massive 



amplitudes were identical to ones appearing in massless amplitudes. An encouraging feature of our 
explicit calculations is that intermediate expressions did not grow in comparison to the size of final 
results. 

Unitarity methods may also be applied to higher-point massive calculations; as an example we 
have presented the contribution of a massive fermion loop to the five-gluon amplitude with identical 
helicities. For general processes one would need to extend the arguments presented in this paper. 
This type of investigation of the cut structure of integral functions for an arbitrary number of 
external legs has already been carried out for the massless case |7]]. 

The unitarity methods discussed in this paper may be combined with other techniques. For 
example, for five- and higher-point amplitudes, factorization provides a strong constraint |L6|, CO. H, 



17]. Supersymmetry [1C, 11] and string motivated [13, 14] ideas can also be incorporated into the 



unitarity approach. One should also be able to devise automated amplitude evaluation programs 



1 32] based on unitarity. 

The same types of methods as those discussed in this paper can be applied to cases where there 
are a variety of masses in the loop, such as for heavy fermion production, and will be discussed 
elsewhere p3f . One should be able to use unitarity as a tool for constructing two- loop amplitudes 



1 24], but to perform four- or higher-point calculations one would first need to construct a table of 



two-loop integrals [34] 
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Figure 2: Examples of labels on integral functions with four-point kinematics. The triangle 1^ in 
(b) is obtained by removing the internal propagator between external legs 1 and 2 in (a). Similarly, 
the bubble I^'^ in (c) is obtained by removing the internal propagator between legs 3 and 4 in 
(b). 



A Integrals 

In this appendix we collect expressions for the integrals [2^, 23] we use in this paper. The 
D-dimensional loop integrals considered in this paper are defined by 

f d D P P ai ■ ■ ■ P a ™ 
I°[P ai ...P Q ™1 =i(-l) n+1 Uir) D/2 -r-rn 1 , (A- 1 ) 

nl i I ^D/n9 ___2\ UP — S^n-l r..\2 _ ^n2\ ' V ' 



where P is a D dimensional momentum and the external momenta, ki, are four-dimensional. Our 
convention is to suppress the dimension label on I n when we are dealing with a D = 4 — 2e integral. 
We also suppress the square brackets when the argument is unity. 

Given a four-point amplitude, we have two conventions for labeling the one-, two- and three- 
point integrals. The first, and more familiar one is to explicitly give the kinematic invariant upon 
which the integral depends. For example, his) is the (4 — 2e)-dimensional bubble that has an 
invariant mass square of s = {k\ + A^) 2 flowing through its external legs. The second convention 
for labeling the lower point integrals is to indicate, by a raised (i) that the internal propagator 
between legs i — 1 (mod 4) and i has been removed. If two indices are present as in then two 
propagators have been removed. This follows the labeling convention of refs. [^3| and is summarized 
in fig. i 

For this paper it is sufficient to give explicit expressions for the (4 — 2e)-dimensional one, two, 
three and four-point integrals to O(e ), but in some cases the all orders expressions are easy to 
obtain. 

A.l The D = 4 — 2e scalar integral functions 

The one-point function is, in closed form, 

As expected from simple power counting, this is quadratically divergent and is singular for both 
e = and 1. 
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The bubble with an external kinematic invariant s is 



I 2 (s) = lf A) = / 2 (0) + 2 + xlog 



x — 1 

X + 1 



+ 0(e) 



(A.3) 



where x = a/1 — 4m 2 /s. For external massless kinematics, this may be written in closed form as 



as 



/ 2 (0) = m 



- 2e r(i + e) 



(A.4) 



The scalar triangle integral is 



h(s)=L 



(2) 



1, 2 f x + 1 
77" lo g 



2,s 



x — 1 



while the box with a uniform internal mass is 



U(s,t) = — 
st 



/ urn? m 2 \ / urn,- m~ 
V st ' s > V st 7 t 



i i 
um m 



+ 0(e) 



where u 



-s — t, 



fr(x,y) = 



x+ — x_ 



Li 2 



Li 2 



x--y 



+ Li 5 



.X i 



2?+ - y 



+ Li 2 



y - x_ 



(A.5) 



(A.6) 



(A.7) 



with x± = i(i± VI -4X), y= 1(1 + V1-4Y), 

and the dilogarithm [35] is Li 2 (x) = — Jq 1 (it ln(l — xt)/t. 

For the case of a single massive external leg, with an invariant mass-square M 2 , we have 



h(s,t,M 2 



1 

st 



2 2 



st 



+ 



2 2 

wm m 



st t 



um 2 m 2 \ 
"^T'M 2 / 



+ (Ai 



where u = M 2 — s — t. Through O(e ), the box integrals contained in Kjp in eq. (|6.4D may be 
obtained from here and from eqs. (|A.12j ) and ( |A.25 ); for example, ij = Zi(s23) s 34! s 5i)- 

The scalar pentagon integral appearing in five-point amplitudes may be obtained from four- 



point integrals using the recursion formula [22 



where 



i=l 



°i ~ S ij 1 > 
i=i 



(A.9) 



(A.10) 



(i) 

and the I\ are boxes with a uniform internal mass and a single massive external leg. The matrix 
Sij is 

Sij =m 2 - ^p 2 j , with pa = , and pij = pji = h + k i+ i H h kj-i for i < j . (A. 11) 

Similar formulas exist for arbitrary massive or massless kinematics [23]. 
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A. 2 The higher dimension integrals 



The higher dimension integrals are defined in eq. ( |A.1| ) with D set to the appropriate value. The 
relationship of these integrals to the usual four-dimensional ones has been extensively discussed 
in ref. [23]. For n < 6 they can be written in terms of the (4 — 2e)-dimensional integrals via the 
integral recursion relations 



r D=6-2e 



1 



-2e 



(n - 5 + 2e)c 
1 

(n - 7 + 2e)c 



21. 



D=4-2e 



i=l 



(i),D=A-2e 
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D=6-2e 



n 



1=1 



(i),D=6-2e 



(A.12) 



where c% is defined in eq. ( A.lOj ) (with the five-point kinematics replaced by n-point kinematics) 
and c = Ya=\ Ci- 

Higher dimension integrals arise naturally when performing the cut calculations described in 
this paper. In a typical cut calculation discussed in the text we obtain integrals of the form 



i(-ir +i (4^) 



\2-e 



(fp cT 2e /i 



(2vr) 4 (2vr)- 2e (p2 



m 2 ) . . . ((p - Yl 



n-l 



fi 2 — m 2 ) 



(A.13) 

where we have explicitly broken the (4 — 2e)-dimensional momentum into a four-dimensional part, 
p, and a (— 2e)-dimensional part, [i. 

Since we are using the Minkowski metric of negative signature and the fractional vector has 
only space-like components, we define fi ■ fi = — fi 2 . Note that dimensional regularization can be 
thought of as altering the mass appearing in the loop propagators. This mass is integrated over 
when we perform the fractional dimensional vector integration: J d~ 2e fj,. 

In general, odd powers of fi a cancel from one-loop integrals. The \x a integration, of eq. ( |A.13 ) 
is formally in a sub-space that does not overlap with any other of the momenta associated with 
the loop. For this reason any contribution to the numerator that is odd in the vector fi a will 
not contribute to the integral. Accordingly, we shall only need to consider the cases where the 
integrand depends on \i a only through fi 2 . Note that we do not need to consider the more general 
cases of \A a pr etc., because in a full amplitude all of the (— 2e)-dimension Lorentz indices must be 
contracted against each other. 

Consider an integral of the form 

' r ~ 2 ' P (^ 2 )7(P Q ,^ 2 ), (A.14) 



(27r) 4 - 2e 

where p and jj, are the four- and (— 2e)-dimensional components of P. Breaking up the measure 
into the two types of components we have 



d 



-2c 



d 4 p 
{2tt) 4 J (2vr) 



(» 2 ) r f(p a ,^ 



-2c 



dfi 2 



{2ir) 2r J dtt- 
J dQ,2r-l-2e 

-e(l-e)(2 



l-2e 



2(2vr)- 2e 
d A p f d 2r - 2 " 



OA 



-l-e+r 



(2vr) 4 
e) • • • (r - 1 



I* 



(27T) 



2r-2e 



(47T) 



d 4+ 



2r-2e 



P 



(2VT) 



4+2r-2e 



3^7 f( P a ,^), 



(A.15) 
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where we used the independence of the integrand on the orientation of the vector n a and the 
formula for solid angles in general dimensions 



/ dn n = f 



In particular, using the definition (A.l), we find 



I, 



£>=4-2er..2 



-el: 



D=6-2e 



and 



r_D=4-2er..4 



[/!*] = -6(1 - 6)4 



D=8-2e 



(A.16) 



(A.17) 



Note that although the loop momentum is shifted to higher dimension, the external momenta 
remain in four dimensions. 



An analogous method leads to the result, 



fftf,i?) = (-4*V (r + 1)! / * Jt /faV) 



(2vr) 4 (2tt)- 2 ^ V V v rj ( 27r ) 4 +2r-2e 

which is useful when evaluating Feynman parameter integrals. 



(A.18) 



A. 3 Feynman parameter reduction 

As discussed in the text, tensor integrals may be evaluated using Feynman parameters. We make 
use of the following formulae in the text. 

After Feynman parameterization and the normal shift of momentum we make use of the fol- 
lowing, 

d^q d~ 2e fi 



7 Q 1 . . . n a 2r + l 



d A q d- 2e fi 

(2vr) 4 (2vr)- 2e {q 2 - fi 2 - S^a^Y 



(2tt) 4 (2tt)- 2 ^ (q 2 - ix 2 - SijUiaj] 
fql 3 1 f d A q (T 2e /x 



0. 



(2tt) 4 (27r)- 2e (q 2 - p? - SijaidjY 

l_ a/3 jD=6-2e 

2 (4vr) 2 - e ' n 



(A.19) 
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where Sij is the matrix defined in eq. ( A.llj) with n-point kinematics. 

After integrating out the loop momenta we obtain integrals of the form 



#[/(»*)] = T ( n ~ D / 2 ) I dna k S(l ~ Er«r) 

J 



f{a k ) 



n-D/2 



(A.21) 



where Sij is defined in eq. ( A. 10] ). An integral reduction formula that is quite useful is [23] 

i n 

T D \ a -] - - Vr • T {j) ' D + — I D 



CO 



where 



c _i c i c j 
^ Co 



(A.22) 



(A.23) 



and ci is defined in eq. ( A.10| ) (with the five-point kinematics replaced by n-point kinematics) and 

En 
i=l c »- 
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A. 4 Common integral combinations 



The following combination of integrals often appears in the amplitudes of this paper. 

J n = I n [m 2 + fj 2 } = I n [m 2 ] + I n [/x 2 ] 

2 T D=4-2e 



and 



K n = I n [(m z + u 2 ) 2 ] = m 4 / n [l] + 2m 2 I n [^} + I n [ U 4 ] 



D=6-2e 
n 



[1] 



4 T 

m L 



4 T D=A-2e 



[l]-2m 2 e/ n c = 6 - 2e [l]-6(l- e )^ 



D=8-2t 



[1] 



(A.24) 



(A.25) 



Using eq. ( A, 12 ) we can remove the explicit e and m 2 dependence from these integral combina- 
tions. Specifically, we find 



J2{s) = \h{s) + \h 

Hs) = \l 2 {s)-L 



D=6-2e 



T D=6-2e 
L 2 

00, 

t 



(A.26) 



Au 



2u 



2u 



hit) 



r D=6-2e 



and 



r _ 3 T D=8-2e 1 s T D=6~2( 
= 7J4 + -T- J3 

4 4 4u d 



4u 



2n 



4u' 



The forms in eqs. ( A.26|) and (A. 27) are more useful for identifying the quadratic divergences. 



(A.27) 



B Fixing remaining ambiguities 

In this appendix we show that the two-step procedure, given in section Q.2, uniquely fixes all cut- 
free functions in the four-gluon amplitudes with massive loops. Here, we fix the coefficients of 
these functions using the known ultraviolet divergences of amplitudes. In general, infrared mass 
singularities also provide a means for fixing such coefficients. 

First we establish the set of cut free functions that may appear in the amplitudes and then 
we explain why the two steps, in section |].2, fix the coefficients of all cut-free functions. For null 
momenta, k 2 — ► 0, the following three functions have no cuts in kinematic variables 

h, h(k 2 ), and / 2 D=6 - 2e (fc 2 ), (B.l) 



where the integral functions I\ and I2 are scalar tadpole and bubble functions in eqs. ( A.2[) , (|A.4|) 



and ( A.12| ). In the massless case (m = 0) these three cut-free integrals identically vanish in 
dimensional regularization and are irrelevant. Observe also that, for k 2 7^ 0, the two-point functions 
do have a cut in k 2 , which means that only coefficients of the tadpoles are cut-ambiguous. Thus, 
there are fewer cut-ambiguous functions to consider for massive external legs. 

First we show that there are no combinations of other integrals which are cut free. To do this we 
must establish the set of functions that can appear in the amplitude. Using a Passarino-Veltman 
(PV) type reduction [36, 21], one may re-express all tensor integrals from the Feynman diagrams 



as linear combinations of D = 4 — 2e scalar integrals. Accordingly, we only need to consider scalar 
integral functions. 
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In general the coefficients of integral functions appearing in amplitudes may depend on e. An e 
in a coefficient of a divergent integral may lead to a rational function at (D(e ), as in eq. 

HI)- if 



the integral function combinations are cut-free, this would complicate our analysis since one would 
need to track the e-dependence. For example, if a coefficient were e dependent one might obtain 
cut-free functions of the form 

C(e) J 2 (0) = (a + be + • • •) m" 2£ Iil±£l = a m~ 2e Eil±f) + 6 + 0(e) . (B.2) 

Whilst the coefficient a could be fixed from knowledge of the ultraviolet divergence (or from m — > 
properties) of the amplitude, without a detailed understanding of the coefficient C(e) we would 
not be able to fix the rational function, b. It is thus convenient to eliminate all e-dependence from 
the coefficients. To do so, we use an appropriate regularization scheme and a slightly modified 
reduction algorithm. 

In addition to altering the number of loop momentum dimensions, conventional dimensional 
regularization also changes the number of physical states; this introduces e-dependence at the level 



of the Feynman diagrams. Using the four-dimensional helicity scheme p8| , 13, [T^], however, the 
number of physical states is the same as in four-dimensions, so no e-dependence is introduced in 
the diagrams. 

We also use a systematic procedure for reducing tensor integrals which produces no e-dependence 
in the scalar integral coefficients^. This modified PV reduction is presented in the first appendix 
of ref. |jl7f l; all e-dependence from the integral reduction is absorbed into higher-dimension scalar 
integral functions. Using the representation for J n and K n in eqs. flAjD and ([AT27D , the four- 



point amplitudes given in this paper may be explicitly expressed in this form. (At any point, the 
amplitudes may be re-expressed solely in terms of the more conventional D = 4 — 2e integrals via 
the integral recursion formulas in eq. ( A. 120 .) 



We need to understand which linear combinations of integrals in a gedanken calculation of an 
amplitude can produce a cut-free function. For four-point amplitudes, this can be found from 
a direct inspection of the integral functions in appendix |A|. (This analysis is similar to the one 
performed for massless amplitudes in ref. |j]].) 

Putting aside the higher-dimensional integrals, we first consider the D = 4 — 2e integrals in the 
first part of appendix |A[ The coefficients of those integrals with kinematic dependence are fixed by 
the cuts because each of the functions contains a logarithm or dilogarithm unique to it: the bubbles 
contain single logarithms, the triangles squares of logarithms, and the box dilogarithms that tangle 
the s and t channels. Consequently no linear combination of these integrals can be formed which 
is cut free. 

To show that the higher dimension integrals with nontrivial kinematics can be distinguished 
with the all orders in e cuts, we make use of the integral recursion relation eq. ( |A.12 ). This 



relationship between higher dimension integrals and D = 4 — 2e integrals introduces a distinct 



e-dependence to coefficients. For example, by eliminating J 4 6 2e in favor of D = 4 — 2e integrals 



with eq. ( A.12|) , a factor of (— 1 + 2e) 1 is introduced into the coefficient of the D = 4 — 2e box. 



In reducing a D = 8 — 2e dimensional box to D = 4 — 2e dimensional ones, an additional factor 
of (—3 + 2c)" 1 is obtained. Since the cuts are computed to all orders in e, the cuts detect the 
full analytic structure of these coefficients. Similarly, to distinguish between I®~ 6 ~ 2e ; I^~ e ~ 2e and 

6 In the Passarino-Veltman paper ]36|, this e-dependence is not explicit. Instead, it appears as constant terms not 
proportional to loop integrals. For example the 1/4 in the reduction function, C24, on pl99 of this paper arises from 
an e in a coefficient multiplying an ultraviolet divergent integral. 
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J^ )=4_2e , we note that after applying the integral recursion formula (A. 12), we obtain D = 4 — 2e 
scalar bubbles, proportional to (—2 + 2e) _1 , (—3 + 2e)~ l and 1, respectively. Given an amplitude 
expressed as a linear combination of D = 4 — 2e and higher dimension integrals, with coefficients 
having no e-dependence, the coefficients of the higher dimension integrals with nontrivial kinematics 
are therefore uniquely specified by the cuts. 

Thus, the only functions with no cuts which may appear in an amplitude are the ones listed in 
eq. ( [B.lD , so any four-point gauge theory amplitude will be of the form 

A A = Cut-constructed part + lim (di h + d 2 I 2 =6 ~ 2e (k 2 ) + d 3 I 2 (k 2 )) . (B.3) 

The coefficients di may contain a 1/k 2 so it is convenient to to write the cut ambiguous parts as a 
limit. A single power of 1/k 2 is the worst that one encounters in a one-loop gauge theory calculation. 
(For a mass in the loop the k 2 — > limit is smooth and no infrared divergences develop.) 

Expanding the integral 1®~ 6_2e [l], defined in eq. ( A.l ), with respect to the external momentum 



we have 

I°=^{k 2 ) = m 2-2,r(l + £l + l fc 2 m -2e£(l + e) + 4 

<e~l) 6 6 (B4) 

= /i + ^ 2 / 2 (0) + O(A; 4 ), 

so I^ )_6_2e (0) contributes only as a linear combinations of 1% and ^(O). Thus, there are only two 
independent cut-free integral functions which may enter an amplitude with a uniform mass in the 
loop. Hence, the amplitude must have the form 

A± = Cut-constructed part + &\ I\ + d 2 ^(0) . (B.5) 

Since there are only two cut-free functions, I\ and /2(0), and two types of ultraviolet divergences, 
quadratic and logarithmic, we may fix the coefficients d\ and d 2 by making the amplitude have the 
correct total ultraviolet divergence. This explains the two step procedure of section ^. 

Following the discussion for the massless case in ref. [jjj one may extend the above arguments 
to an arbitrary number of legs. (One difference is that when working to all orders in e, scalar 
pentagon integrals can no longer be eliminated from the set of integrals in terms of which the 
answer is expressed.) 



C Rules for fi a and \ji 

The identification of dimensional regularization with the introduction of an integrated mass-like 
vector, /i Q , necessitates some discussion. To begin, we have a (4 — 2e)-dimensional vector, Q a , 
which we can express as a sum of a four dimensional piece q a and a remainder /i a , Q a = q a + /x a . 
It follows that Q = d + Here we discuss the properties of ft. 

Firstly, some definitions. We abide by the usual conventions for the Dirac algebra, 

{ 7 « /} = 2r] a P , 7 Qt 7° = 7 V . (C.l) 

In this equation, a and (3 are (4 — 2e)-dimensional Lorentz indices and the metric is rj a ^ = 
diag(+, — , — , .). It follows that, 

2Q-Q = = {$M} + {^}+2{^} =2q- q + 2w + 4q ■ = 2(q 2 - /j 2 ) . (C.2) 
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The cross-term vanishes because we have formally chosen fi a to be in a sub-space orthogonal]] to 
the four-dimensional subspace containing q a . Without this choice, terms of the form fi ■ k would 
exist, significantly complicating our calculations. The minus sign in going from fj, ■ fj, to /U 2 is from 
the metric. For the familiar four-dimensional vector, as is conventional, we write q ■ q = q 2 . 

When the (4— 2e)-dimensional vector is null, we find that the four-dimensional vector effectively 
becomes 'massive': 

Q-Q = 0^q 2 = fi 2 , (C.3) 

where we denote a (4 — 2e)-dimensional momentum by a capital letter and a four-dimensional 
momentum by a lower-case letter. The quantity fj? is always integrated over. 

Since the metric is diagonal, d freely anti- commutes with four-dimensional j a matrices, 



= 0. (C.4) 

This has an interesting consequence with respect to manipulations with 75. We use the conventions 
of 't Hooft and Veltman [29, 27|, and adopt the arbitrary dimension definition, 75 = Z7°7 1 7 2 7 3 . As 
usual, 75 anti-commutes with all the four-dimensional Dirac matrices: 

,75} = 0. (C.5) 



With respect to jl, however, 75 freely commutes, 



0, (C.6) 



since 75 is defined in this prescription as a product of an even number of four-dimensional Dirac ma- 
trices. The (— 2e)-dimensional component, d, commutes freely with the helicity projection operator, 
u± = i(l ±75), 

UJ ±^ = ^ LJ T^ but, w± d = duj± . (C7) 

In our dimensionally regularized cut calculations the momentum of cut fermion lines is to be 
considered (4 — 2e) dimensional. We will borrow the bra and ket symbols to represent the spinors, 
but as helicity is not a good quantum number we shall not label them with a helicity. (Since we 
always sum over all states across the cut, there is anyway no need to define a (4 — 2e) helicity 
notion.) These spinors obey the conventional Dirac equation, 

(Q\Q = (Q\m, {Q\d=(Q\{m-if). (C.8) 

To handle trailing fermions we will use the same expression but in such a case the momenta will 
be explicitly negative, d\ — Q) = —{m + d)\ — Q). We adhere to the convention that the argument 
momentum is the momentum flowing in the direction of the fermion arrow. 

In sewing the (4 — 2e)-spinors together (across the cut) we implicitly sum over the two spin 
degrees of freedom, that is 

\Q)(Q\=Q + m = d + i/L + m, and | - Q){-Q\ = -$ + m = -<| - ^ + m . (C.9) 

As mentioned in the text, this notation glosses over the distinction of spinors and anti-spinors and 
can introduce an overall minus sign that needs to be put back by hand. 



7 This choice does not correspond to the one made in the dimensional reduction [Eq] regularization scheme. 
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Having kept a four-dimensional definition for 75, we can apply spinor-helicity methods (8| to 
compress tree amplitudes. However, some care is required in applying the identity ( |2,4[) . Whereas 
this is the conventional rule for substituting for a slashed polarization vector, it is only a valid iden- 
tity when the Dirac algebra on at least one side is four-dimensional. For example, for (Q\^d^d . . ., 
before applying the rule (|2.4j ) we must first anti-commute each of the d to one end of the spinor 
line. In this way we pick-up factors of — dd = /u 2 , just as we would combine conventional mass 
contributions from the fermion propagator. Once we have two neighboring polarization vectors, 
^ ± (k;r)^ ± (p; q), we are free to apply the rule Q2.4|). 
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